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Abstract 

A quantum-dynamical analysis of phonon-driven exciton dissociation at polymer hetero- 
junctions is presented, using a hierarchical electron-phonon model parameterized for three 
electronic states and 24 vibrational modes. Two interfering decay pathways are identi- 
fied: a direct charge separation, and an indirect pathway via an intermediate bridge state. 
Both pathways depend critically on the dynamical interplay of high-frequency C=C stretch 
modes and low-frequency ring-torsional modes. The ultrafast, highly non-equilibrium dy- 
namics is consistent with time-resolved spectroscopic observations. 
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Semiconducting 7r-conjugated polymers are promising low-cost and flexible materials for 
electronic devices such as organic light-emitting diodes (OLEDs) and solar cells [1-3] . Many 
of the unique properties of these materials stem from the combination of the molecular 
structure of their building blocks and the extended nature of the vr-system. Thus, excitonic 
states, i.e., electron-hole quasi-particle states that extend over the polymer strands, play 
a central role for the optoelectronic properties [3]. Exciton formation can result from 
polaron recombination, e.g., in OLEDs, or else from the primary photoexcitation step, 
e.g., in photovoltaic systems. 

A key feature of organic semiconductors is the strong coupling between the elementary 
electronic processes and the material's phonon modes [3, 4] - much stronger than in inor- 
ganic semiconductors. This is underscored by experimental observations of sub-picosecond 
scale coherent nuclear motions following photoexcitation [5] . Further, the ultrafast nature 
of the electronic decay processes [6-8] suggests the presence of coherent vibronic coupling 
mechanisms, possibly determined by conical intersection topologies - very similar to the 
photophysics of polyatomic molecules and to Jahn- Teller effects in sohds [9, 10]. These 
observations call for a microscopic-level, quantum-dynamical interpretation of the basic 
charge and energy transfer processes. 

In this Letter, we focus on the role of electron-phonon coupling in the charge transfer pro- 
cesses at bulk polymer heteroj unctions [11]; these provide conditions for an extremely effi- 
cient charge separation at the interface between different phase-segregated polymers. The 
primary excitation is a photogenerated exciton stabilized by the electron-hole Coulombic 
interaction (with a typical binding energy of eg ~ 0.5 eV [2]). Due to the highly folded 
interfacial area in bulk heterojunctions, the exciton has a high probability of reaching the 
interface within the diffusion length (typically ~ 20 nm). The exciton decay towards a 
charge-separated state ("exciplex") is largely determined by molecular-level electronic in- 
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teractions at the interface. To a first approximation, the efficiency of charge generation 
depends on the ratio between eb and the band offset between the two polymer species 
[11,12]. Recent time-resolved photoluminescence studies have shown that the exciton de- 
cay - and possible exciton regeneration - fall into a (sub) picosecond regime [7,8]. While 
the elementary processes at the polymer interface are crucial for the device function and 
optimization, these processes are not well characterized as yet on the microscopic side. 
Against this background, the purpose of this Letter is to present a realistic quantum- 
dynamical model involving several electronic states and an explicit representation of the 
phonon distribution. 

We specifically consider a hetero junction composed of the poly[9,9-dioctylfiuorene- 
co-bis-N,N-(4-butyl-phenyl)-bis-N,N-phenyl-l,4-phenylenediamine] (TFB) and poly[9,9- 
dioctylfiuorene-co-benzothiadiazole] (F8BT) polymer components [7, 12-14]. At the 
TFB:F8BT hetero junction, is found to be similar to the band offset, so that both exciton 
decay and regeneration phenomena are expected to occur. Semi-empirical calculations were 
used to identify the relevant electronic states and parameterize the vibronic coupling Hamil- 
tonian used below (see Eq. (1)) for 24 explicit phonon modes covering the high-frequency 
(C=C stretch) and low-frequency (ring-torsional) phonon branches [13,15,16]. The exciton 
state (XT), which is the lowest excited state with significant oscillator strength, mainly 
has electron-hole density on the F8BT moiety, while the lowest-lying interfacial charge 
transfer state (CT) corresponds to a TFB"'"-F8BT~ charge separated state. (However, the 
CT state is not purely polaronic since it does carry oscillator strength to the ground state 
through mixing with excitonic configurations on the F8BT moiety). Among the manifold 
of remaining states, one is found to exhibit a strong coupling to the XT state [13]; this 
state will be considered as an intermediate state (IS) in the following, see Fig. 1. 

To interpret the nonadiabatic decay of the exciton state, we apply a three-state vibronic 
coupling model in conjunction with highly accurate wavepacket propagation techniques. 
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i.e., the multiconfiguration time-dependent Hartree (MCTDH) method [17]. In particular, 
we use a hnear vibronic coupUng (LVC) model [9] which is appropriate in view of the 
moderate nuclear displacements during the process. The A'"-mode LVC Hamiltonian (here, 
N—2A) reads as follows, using mass and frequency weighted coordinates and atomic units 
{h = 1), 



N N 
i=l i=l 



[ k!^x, )S^x^ X^'^x, \ 



Xi X{ Xi 

I ^ (13) . (23) (3) 
\ "^i '^i '^i 



+ c 



(1) 



Here, 1 and C denote the unit matrix and a coordinate independent constant matrix, re- 
spectively; the uji, Pi — —id/dxi, and Xi are the frequencies, momenta, and displacements 
along the vibrational normal modes. The diagonal and off-diagonal potential terms corre- 
spond to the diabatic potentials and couplings. The model Eq. (1) allows for the presence 
of conical intersections between pairs of electronic states at nuclear configurations where 
the diabatic couphng vanishes and the adiabatic states become degenerate [9] . 

We further apply a recently developed effective-mode representation of the LVC model 
[18], which separates the Hamiltonian into effective-mode vs. residual- mode parts. 



H — HeS + H 



(2) 



By an orthogonal coordinate transformation, X — Tx, a subset of effective modes is 
generated which subsume all vibronic coupling terms of an electronic A^^-state system into 
\Ns{Ns + 1) effective modes [18]. This construction, which was previously applied to two- 
state systems (three effective modes) [16, 18, 19], is here extended for the first time to three 
electronic states (six effective modes) , 
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where the parameters and modes X are defined by extension of the two-state construction 
described in Refs. [16,18]. Two of the six effective modes (i.e., Xi and X2) are chosen as 
topology-adapted modes which span the branching plane [20] for a given pair of electronic 
states (here, states 1 and 2), i.e., the plane along which the degeneracy at the conical 
intersection is hfted, see Fig. 2. Several of the couphng constants in Eq. (3), i.e., the Af'^'' 

(12) (13) 

(i > 2), Di (i > 3), Ki [i = 4, 5, 6), (^ = 5, 6), and Ag , are zero by construction. 

The residual modes, X^, i — 7, . . . N, do not couple directly to the electronic subsystem, but 
couple bilinearly to the effective modes. By analogy with our recent two-state analysis [16], 
i?res can be transformed to a band-diagonal structure, yielding the nth-order Hamiltonian 

n 

H^^^ = H,s + Yl ^ril (4) 

1=1 

where each /th-order residual term comprises 6 modes, 

61+6 ^ 61+6 i-l 

i=6l+l 1=61+1 j=i-6 

For 6 + 6n — N, the nth-order Hamiltonian Jf^"^ is equivalent to the original Hamiltonian 
Eq. (1). 

We refer to Eqs. (2)-(3) in conjunction with Eqs. (4)-(5) as a hierarchical electron-phonon 
(HEP) model [16, 19] . Successive orders of the HEP Hamiltonian can be shown to con- 
serve the moments (cumulants) of the original Hamiltonian up to the (2n-|-3)rd order [16]. 
Truncation of the HEP Hamiltonian creates a series of reduced dimensionality models, 
associated with approximate nth-order propagators which accurately reproduce the dy- 
namics of the overall system up to a certain time. While the lowest-order description - 
Hcs only - can be appropriate, e.g., for a rapid passage through a conical intersection 
region [18], the present system will be shown to necessitate a higher-order analysis. 

The effective-mode construction for the TFB:F8BT heteroj unction shows that Hcs is es- 
sentially determined by high-frequency C=C stretch modes while J?res determined by 
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low-frequency ring torsional modes. This is analogous to our previous analysis of a two- 
state XT-CT model [16]. Importantly, the low-frequency modes will be shown to play a 
key role in generating the observed ultrafast exciton decay. 

Fig. 2 shows the diabatic XT, CT, and IS potential energy surfaces (PES), for a cut through 
the XT-CT branching plane defined in terms of the effective modes {Xi, X2)xt,ct, and an 
alternative representation involving the XT-IS branching plane spanned by another set of 
effective modes (X^, X2)xt,is- As can be inferred from the figure, the exciton dissociation 
process is determined by a landscape of multiple intersecting surfaces, and the branching 
plane topology is a key factor in the nonadiabatic dynamics. Of particular interest for the 
present analysis is the role of the intermediate state (IS) that lies above the exciton state 
but could be accessible due to its strong diabatic coupling to the XT state, and due to the 
significant extension of the wavepacket. 

Fig. 3a shows the population evolution obtained from MCTDH simulations for the overall 
3-state, 24-mode system according to Eq. (1). An ultrafast initial XT CT decay takes 
place (~ 50% at 150 fs), followed by an oscillatory behavior of the XT and CT populations 
which is found to persist over the 3 picosecond observation interval. These observations are 
in qualitative agreement with the experimentally observed exciton decay and regeneration 
[7] , even though a steady overall increase of the CT population is eventually observed in 
our simulation. The IS population experiences an initial increase, but does not rise above 
an average of ~ 10% beyond 500 fs. Yet, one cannot exclude the possibihty that this state 
acts as a "bridge" facihtating the XT-CT transfer. 

To investigate the possibility of an exciton decay pathway via the IS state, XT — > IS — > 
CT, we set the XT-CT diabatic coupling artificially to zero, such that only the indirect 
pathway is permitted. As can be inferred from Fig. 3b, a subpicosecond exciton decay 
and concomitant rise of the CT state population is also observed in this case, even though 
the process is somewhat slower than in the fully coupled 3-state dynamics. Again, the IS 
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population tends towards an average of ~ 10%, i.e., the IS state facilitates the XT-CT 
transition but does not act as an additional acceptor state. Overall, the dynamics of Fig 
2a is therefore the result of a superposition of the direct (XT —>■ CT) and indirect (XT —>■ 
IS — > CT) pathways. Given that the XT— > IS transfer is endothermic (with a ~ 0.2 eV 
barrier) in a conventional kinetic picture, the efficiency of the indirect pathway is entirely 
a consequence of the quantum dynamical character of the process. 

The remaining part of the analysis addresses the HEP construction of Eqs. (2)- (5), with 
the twofold objective of identifying a reduced-dimensionality picture of the dynamics, 
and understanding the role of the different types of phonon modes, i.e., high-frequency 
C=C stretch modes (constituting Hes) vs. low-frequency ring-torsional modes (constitut- 
ing -f^res)- As shown in Fig. 4, the high-frequency modes by themselves, at the level of 
the 6-mode H^^^ approximation, cannot reproduce the exciton decay (even though they 
correctly reproduce the shortest-time dynamics, on the order of 50 fs). At this level of ap- 
proximation, quantum phase coherence effects are substantially overestimated, leading to 
a weak, periodic population transfer. By contrast, inclusion of the low-frequency modes, at 
the level of the 12-mode H^^^ approximation, results in a qualitatively correct description 
of the dynamics. Indeed, the low-frequency modes induce vibrational energy redistribu- 
tion and dephasing effects that eventually lead to an irreversible exciton decay. Finally, the 
H^'^\ 18-modc approximation is essentially identical to the exact, 24-modc result on a ~ 1 
ps time scale. A similar interpretation holds for the indirect pathway in the absence of XT- 
CT coupling. All of these results would carry over to an initial parameterization involving 
a very large number N of phonon modes, e.g., thousands of modes. While dynamical cal- 
culations according to Eq. (1) would not be feasible for such a "macrosystcm" , the HEP 
model gives a correct reduced-dimensionality description of the ultrafast dynamics. 

In summary, we present the results of fully quantum dynamical calculations of a model 
polymer hetero junction interface. In particular, we highlight the role that intermediate 
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states may play in the overall XT— >CT conversion. While these results go significantly be- 
yond our earlier interpretation for a two-state model [16], the key role of the low-frequency 
modes is confirmed, pointing to the generic nature of the observed dynamical pattern for a 
two-band phonon distribution (see also Refs. [4,21] which address the spectroscopic signa- 
ture of both types of phonon modes in phenylene based polymers) . Our analysis highlights 
the coherent, quantum-dynamical character of the process, and the importance of the mul- 
tidimensional topology. Thus, the branching plane of Fig. 2 replaces the conventional Mar- 
cus parabola picture [22] , and the HEP model provides the correct reduced-dimensionahty 
dynamics on short time scales. Extensions of the present analysis concern (i) the inclusion 
of temperature efi'ects, and (ii) the improvement of the present parameterization, along 
with the generalization to multiple electronic states that potentially appear in the role of 
bridge states. The positions of such states relative to the XT and CT states and their elec- 
tronic properties are indeed quite sensitive to the local morphology and packing geometries 
of the respective polymers at the interface as indicated by recent time-dependent density 
functional theory [14] and semi-empirical investigations [23]. Consequently, the extent to 
which such intermediate states will play a role will depend on the local morphology of the 
interface and the actual physical system will be an ensemble of individual cases. Using 
this local, molecular-level information, the HEP approach should contribute to a realistic 
modeling of the highly nonequilibrium exciton decay processes in extended 7r-conjugated 
systems. 

We thank Audrey Pereverzev, Etienne Gindensperger, and Lorenz Cederbaum for con- 
structive discussions. This work was supported by the ANR-05-NANO-051-02 and ANR- 

NT05-3-42315 projects, by NSF grant CHE-0345324, and by the Robert Welch Foundation. 
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Figure legends 



Fig. 1. (Color online) Schematic illustration of (a) the TFB:F8BT donor-acceptor hetero- 
junction and (b) the photoexcitation process and potential crossings. The blue, red, and 
purple lines indicate the XT, CT, and IS diabatic potentials, respectively. (Dotted lines 
indicate the corresponding adiabatic potentials.) The red and purple arrows indicate the 
direct XT — > CT vs. indirect XT ^ IS — > CT exciton dissociation pathways. 

Fig. 2. (Color online) Branching-plane projections of the coupled diabatic PESs: (a) XT, 
CT and IS PESs projected onto the XT-CT branching plane, and (b) XT and IS PESs 
projected onto the XT-IS branching plane; here, the wavepacket width and trajectory of 
the wavepacket center (blue arrow) are also indicated. (These branching space projections 
are associated with two different sets of effective modes X.) The white, blue, red and 
purple circles indicate the conical intersection, the minima on the XT, CT and IS states, 
respectively. The dotted lines indicate the avoided-crossing seam lines. 

Fig. 3. (Color online) Time-evolving state populations for (a) the full 3-state 24-mode 

wavepacket propagation, and (b) a complementary 3-state 24-mode calculation where the 
indirect (XT ^ IS ^ CT) pathway was selected by setting the XT-CT diabatic coupling 
to zero. 

Fig. 4. (Color online) CT populations for the truncated HEP Hamiltonian H^"''' of Eq. 
(4), for n = (6 modes), n — 1 (12 modes), and n — 2 (18 modes), as compared with the 
exact 24-mode result. 
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